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ABSTRACT 


><ru  *£«"'  <*«»/•  es, 

applied  a Nearly  Free  Electron  Model  (NFEM)  combined 


with  a local  pseudopotential  approximation  to  the  analysis  of  the  high 


Within  this 


obtain* qualitative  agreement  with  the 


experimental  observations  of  three  anisotropic  effects:  self. 


electromigration.  Hall  effect,  and  diffusion  thermoelectric  power, 


The  electromigration  of  a weak  scattering  impurity  in  a Bloch 


electron  gas  is  reconsidered  while  taking  into  account  the  conduction 


electron  screening  in  the  Hartree  approximation.  Our  formal  expression 


for  the  driving  force  depends  explicitly  on  the  dielectric  matrix,  and 


reduces  to  previously  derived  formulas  in  the  appropriate  limits.  An 


approximate  evaluation  of  the  formulas  yields  values  for  Z 


Anisotroulc  relaxation  times 


from  the  observed  phonon  dispersion  curves,  tend  to  favor  mass  transport 


in  the  basal  plane 


The  Hall  effect  is  analyzed  in  terms  of  a model  Fermi  surface 
based  upon  the  standard  2-01 V (Orthogonal! zed  Plane  Wave)  approximation, 
The  effect  of  k dependence  of  the  relaxation  times  is  found  to  be 


negligible  in  the  calculation  of  R(J  , but  is  probably  important  for 
R,  . As  temperature  increases,  the  changing  c/a  ratio  and  decreasing 


band  gaps,  as  affected  by  the  Debye-Waller  factor,  combine  to  reduce 


the  values  of  the  Hall  constants.  Our  predictions  for  R„  and  R 


appear  to  agree  quite  satisfactorily  with  values  extrapolated  from 


measurements  at  lower  temperatures 


Tii 


We  simplify  the  Fermi  surface  model  further  in  our  treatment 
of  thermoelectric  power.  As  in  the  noble  metals,  the  shape  of  the  foxm 
factor  results  in  a positive  contribution.  Anisotropy  in  our  calculated 
values  of  S and  s is  due  primarily  to  distortion  of  the  Fermi  surface. 
Our  predicted  values  appear  reasonable  in  light  of  available  experimental 
data  at  lower  temperatures. 


INTRODUCTION 


A.  Historical  Revlev 

1.  General  Transport  Theory 

The  classical  theory  of  electron  transport  in  metals  was 
1 2 

developed  by  Drude  and  Lorentz  in  terms  of  a relaxation  time  that 
represented  the  mean  time  between  collisions.  Sommerfeld^  reformulated 
the  problem  in  a quantum  mechanical  framework.  He  used  Fermi-Dirac 
statistics  for  the  electrons,  but  he  retained  the  concept  of  a simple 
relaxation  time  to  describe  the  effect  of  collision  processes. 

More  recent  transport  theory  has  been  based  upon  the  linearized 
Boltzmann  equation  (Eq.  1.1),  which  describes  the  balance  at  steady 
state  between  the  effects  of  the  applied  electric  field  E?  and  conduction 
electron  scattering. 


(1.1) 


Here  g(k)  represents  the  steady  state  deviation  of  the  electron 
distribution  from  the  the  Fermi-Dirac  function,  f°.  W (£,"£)  is 
proportional  to  the  transition  probability  from  state  k to  It*  and  C and 
v are  the  electron  energy  and  velocity.  The  relaxation  time  approximation 
consists  of  taking  the  function  g(K)  to  be  of  the  form 


(1.2) 


where  again  *v(1t)  is  the  electron  velocity.  In  the  case  of  free  electrons ; 
l.e. , for  a spherical  Fermi  surface  and  w($/k}  depending  only  on  |k-k^  , 

It  can  rigorously  be  shown"*  that  a constant  relaxation  time  exists  in 
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I 

I 

I 

I 


the  form 

i - j’[l-cos©(^W(k,k^#  . 

Fermi  Surface 

This  model  satisfactorily  describes  the  resistivity  of  the  monovalent 
metals,  but  does  rather  poorly  with  the  multi valents.  Obviously,  it 
cannot  hope  to  deal  with  anisotropies  in  non-cubic  systems,  the  problem 
to  which  the  present  work  is  directed. 

Several  workers  have  attempted  to  modify  Eq.  (1.3)  in  some 

physically  reasonable  way 11  For  example,  Giaever0  used  Eq.  (1.3)  to 

define  a ^-dependent  "relaxation  time."  Truant  and  Carbotte11  later  used 

the  same  formula  to  calculate  a relaxation  time  for  Zn,  ignoring  tne 

anisotropic  structure  of  the  material.  Such  approximations  met  with 

limited  success,  but  they  can  result  in  considerable  error,  even  in  the 

12 

relatively  simple  noble  metal  systems. 

One  possible  approach  is  to  avoid  the  direct  solution  of 
Eq.  (1.1) , solving  directly  instead  for  the  resistivity  using  a 
variational  procedure^  This  technique  has  been  applied  with  varying 
degrees  of  success  to  many  metals^-1^  However,  as  pointed  out  by 
Robinson  and  Dow,'  variational  ratios  (e.g.  for  the  resistivity)  which 
explicitly  take  anisotropy  into  account  do  not  simplify  sufficiently  to 
appear  advantageous.  Also,  the  presence  of  anisotropy  significantly 
increases  the  ambiguity  and  uncertainty  inherent  in  the  choice  of  tried 
functions. 

A general  method  of  attewking  Eq.  (1.1)  is  to  define  g(k)  in 
terms  of  a vector  mean  free  path  T(k)}^  i.e. , 

g(Jt)  «■  -(g^-)eE»T(Jc)  . (I*1*) 


! 
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I 

I 

I 
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19 

Taylor  discussed  the  solution  of  Eq.  (l.l)  directly  for  the  vector 

t20 
. It  is  also  possible  to  introduce  a modified  "relaxation  time"  by 

the  relation 

X(k)  =Tv($)  (1.5) 

where  is  in  general  a second  rank  tensor.  The  three  independent 

quantities  determined  by  the  Boltzmann  equation  can  be  taken  to  be  the 

diagonal  components  of  7"  referred  to  its  principal  axes.  Chan  and 
21-22 

Huntington  used  this  approach,  combined  with  a Nearly  Free  Electron 

23-24 

Model  (NFEM)  based  upon  the  pseudopotential  concept  to  investigate 

the  conductivity  anisotropy  in  single  crystal  Zn.  The  Chan,  Huntington 

calculation  explained  how  large  elastic  and  Fermi  surface  anisotropies 

o c 

in  Zn  combine  to  result  in  a nearly  isotropic  conductivity  tensor. 

The  success  of  this  calculation  suggested  using  a similar  model,  to 
investigate  other  transport  parameter  anisotropies. 

2.  Driving  Force  for  Electro.migration 

One  of  the  transport  phenomena  sensitive  to  anisotropies  in 
the  metal  system  is  the  driving  force  for  electronigration,  or  mass 
transport  Tn  a metal  carrying  a high  electric  current  density.  The 
force  itself  is  derived  from  the  measured  drift  velocity  through  the 
Nernst-Einstein  relation,  and  is  usually  expressed  in  terms  of  an 
effective  charge  relating  it  to  the  applied  electric  field.  The  total 
force  on  the  diffusing  atom  has  two  parts,  the  electrostatic  force  of 
the  applied  electric  field,  and  the  "electron-wind”  force,  due  to 
collisions  between  the  atom  and  the  current-carrying  electrons.  We 
shall  be  dealing  primarily  with  the  second  contribution,  that  due  to 
the  streaming  charge  carriers. 
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27  28 

The  early  theories  of  Fiks  and  Huntington  and  Grone 

determine  the  wind  force  by  calculating  the  momentum  transfer  by  the 

electrons  per  unit  time  as  they  are  scattered  by  the  atom.  Bosvieux 
29 

and  Friedel  noted  that  the  electron  current  induces  a dipolar-like 

distortion  in  the  screening  cloud  associated  with  the  diffusing  atom, 

resulting  in  a net  electrostatic  force.  Later,  Sorbello^  extended 

this  viewpoint  to  a more  realistic  model  of  a metal  within  the  framework 

of  pseudopotential  theory.  A somewhat  related  transport- theory  based 

screening  calculation  was  performed  by  Das  and  Peierls?1  More  recently, 

32-34 

a number  of  calculations  have  been  performed  using  Kubo's  linear 

35 

response  formalism.  All  of  these  viewpoints  lead  essentially  to  the 
same  conclusions  concerning  the  wind  force  in  the  FEM;  i.e.,  that  the 
force  depends  on  the  specific  resistivity  of  the  diffusing  atom  and  on 
the  number  of  electrons  per  atom. 

In  order  to  treat  anisotropic  effects  realistically,  the  theory 
must  be  extended  to  include  band  structure  effects.  The  first  efforts 
in  this  direction  were  made  by  FiksJ  " and  by  Feit  and  HuntingtcnT 
In  both  of_these  theories,  the  diffusing  atom  was  considered  to  be  in 
an  extended  plane  wave  state,  an  assumption  which  appears  to  be  more 
appropriate  for  a free  particle;  e.g. , a neutron.  In  addition,  neither 
treatment  dealt  with  the  complications  arising  from  screening  in  a Bloch 
electron  gas.  No  quantitative  evaluations  of  band  structure  effects  in 
electromigration  have  been  published. 

We  are  especially  interested  in  the  effective  charge  of  Zn 

39 

since  it  is  one  of  the  materials  in  which  single  crystal  measurements 
have  been  made.  In  this  work,  Routbort  found  that  the  driving  force  was 


9P!W2221Sl 


5 


about  half  as  effective  along  the  hexagonal  axis  as  perpendicular  to 
it.  In  light  of  the  nearly  isotropic  conductivity  of  Zn,  these  results 

were  Initially  puzzling,  which  stimulated  the  work  of  Chan  and 

22 

Huntington  as  well  as  the  present  work. 

3.  Hall  Effect 

The  transport  properties  of  a solid  may  be  profoundly  altered 
by  the  application  of  a magnetic  field.  Analysis  of  these  effects 
requires  inclusion  of  an  additional  term  in  the  Boltzmann  equation  for 
the  distribution  function  g(lt)  (see  PART  III).  In  the  relaxation  time 
approximation,  Jones  and  Zener^  proposed  a solution  for  g(k)  as  an 
expansion  in  powers  of  the  magnetic  field  ”h.  In  the  weak  field  limit, 
the  resulting  expression  for  the  Hall  coefficient  in  cubic  materials  is 
quite  sensitive  to  the  shape  of  the  Fermi  surface,  involving  the  second 
derivatives  of  the  electron  energy  (first  derivatives  of  electron 
velocity)  with  respect  to  wavevector  k.  Ziman  applied  this  model  to 
a calculation  of  the  Hall  effect  in  the  noble  metals.  Using  what  he 
called  an  "eight-cone  model"  to  account  for  distortions  of  the  Fermi 
surface  from  sphericity,  he  obtained  quite  satisfactory  agreement  with 
experiment. 

18  19 

Price  and  Taylor  have  shown  how  to  extend  the  method  of 
Jones  and  Zener  to  anisotropic  materials  in  terms  of  the  vector  mean 
free  path, *X.  In  this  case,  the  Hall  coefficient  involves  the  derivatives 
of  T with  respect  to  "k.  Calculations  using  the  above  model  have  not 
yet  been  carried  out  for  the  hexagonal  metals. 

Experimental  measurements  of  the  Hall  coefficients  in  Zn  and 
Cd  single  crystals  were  made  by  Lare,  et.  al.  in  1964  and  more  recently 


by  Yonemitsu  and  Takahashl.  They  reveal  a large  anlBOtropy  and 


Indicate  the  predominance  of  "hole-like"  carriers,  while  both  the  a 

and  c-axls  values  tend  to  decrease  with  increasing  temperature.  On 

the  other  hand,  the  electromigration  experiments  suggest  "negative 

charge  carriers";  thus,  there  is  no  clear,  simple  correlation  between 

39 

the  electromigration  results  and  the  Hall  constants.  Recent  calculations 

44 

by  Cowley  and  Stringer  using  the  Plane  Faced  Energy  Surface  (PFES) 

4S-46 

method  achieved  rough  agreement  with  experimental  values  at  room 
temperature. 

4.  Thermoelectric  Power 

The  usual  derivation  of  the  expression  for  the  thermoelectric 
47 

power  of  a metal  Is  easily  generalized  to  anisotropic  metals  (6ee  PART 


IV)  and  gives  the  sane  result 


3 — •ZTST3*' 

- J4»' 


(1.6) 


where  6(0  is  the  electrical  conductivity  of  a metal  whose  Fermi  energy 

is  • Comparison  of  the  theory  with  experiment  is  complicated  by  the 

48 

effect  of  "phonon  drag",  which,  however,  disappears  at  nigh  temperatures 

and  will  be- neglected  here.  Taking  e as  a negative  quantity,  one  finds 

from  the  formula  that  a free-electron  metal  should  have  a negative 

thermopower.  As  is  well  known,  however,  the  thermopover  of  Li  as  well 

as  all  three  noble  metals  is  positive  in  sign. 

4l  49 

Ziman  ^ ^ 0 + famT»+aH  frt  MtaAl va  fM  e wav*o^AV  f a 


attempted  to  resolve  this  paradox  for  the  noble 


metals  by  considering  the  effect  of  Fenni  surface  distortion,  but  was 


unsuccessful.  Taylor  later  modified  the  model  for  Cu  by  Including 


additional  band  structure  effects,  but  succeeded  only  in  driving  his 
theoretical  values  farther  away  from  the  experimental  values.  A new 


7 


' 

\ 


50 

line  of  thought  was  suggested  by  Robinson,  and  developed  to  account 

17 

for  the  positive  thermopower  of  Li.  It  is  that  a critical  determinant 

of  transport  properties  in  the  simple  metals  is  the  shape  of  the  form 

factor  (Fourier  transform  of  the  effective  scattering  pseudopotential) 

resulting  from  interference  between  Coulomb  attraction  and  core  repulsion 

51-53 

potentials.  Subsequent  calculations  for  the  noble  metals  have 

quantitatively  verified  Robinson's  qualitative  observations. 

54  55 

Measurements  * for  Zn  and  Cd  also  show  positive  values  of 
the  thermopower,  as  well  as  a large  anisotropy  between  a and  c-axis 
values.  To  date,  however,  no  calculations,  such  as  the  above  described 
work  on  the  noble  metals,  have  been  carried  out  for  anisotropic  materials. 
B.  Introduction  to  the  Problem 

56 

As  pointed  out  recently  by  Hurd,  the  study  of  low  field 
galvanomagnetic  effects  in  anisotropic  metals  is  presently  a somewhat 
neglected  subject.  Much  of  the  existing  data  is  either  quite  old  or, 
in  many  cases,  incomplete  in  terms  of  temperature  dependence.  On  the 
theoretical  side,  complications  arising  from  the  anisotropic  structure 
and  from  election  band  effects  have  largely  discouraged  quantitative 
investigations.  It  is  the  aim  of  the  work  presented  here  to  take  a 
more  comprehensive  look  into  the  transport  properties  of  Zn  than  has 
heretofore  been  attempted. 

The  model  used  here  for  the  Fermi  surface  and  electron 

vavefunctions  is  consistent  with  that  used  by  Chan  and  Huntington  in 

their  investigation  of  the  high  temperature  conductivity.  It  is  a 2-OFW 

4i 

model,  similar  to  that  employed  by  Ziman  for  the  noble  metals,  in 
which  the  effects  of  the  most  important  Brlllouln  Zone  discontinuities 


. . ^ 
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are  accounted  for,  one  pair  at  a time.  While  this  Is  a highly  approximate 

model,  and  perhaps  more  appropriate  for  the  simpler  noble  metals,  we 

feel  It  retains  the  essential  band  structure  effects  in  the  hexagonal 

metals.  In  addition,  we  assume  throughout  the  existence  of  a local 

pseudopotential  form  factor  (see  APPENDIX  A)  determined  semi-empirically 

57 

from  a prescription  given  by  Stark  and  Falicov.  In  light  of  these  and 
other  approximations  invoked  during  the  course  of  the  following 
calculations,  the  results  presented  here  are  not  expected  to  produce 
exact  agreement  with  experiment;  in  fact,  high  temperature  measurements 
of  the  Hall  effect  and  thermoelectric  power  are  incomplete.  Instead, 
we  hope  to  present  a consistent,  semi-quantitative  picture  of  the 
essential  physical  considerations  resulting  in  transport  parameter 
anisotropy. 

PART  II  concerns  the  driving  force  for  electromigration. 

Pirst,  the  general  theory  of  band  structure  effects  in  electromigration 

la  reconsidered  by  combining,  to  a certain  extent,  the  approaches  of 

Feit^®  and  Sorbello^  The  effects  of  band  structure  on  electron 

screening  are  explicitly  taken  into  account  using  the  "dielectric 
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matrix"  formalism.  * Several  simplifying  approximations  are  made 
and  model  calculations  are  carried  out  for  Zn.  In  PART  III,  the  two 
independent  Hall  constants  of  Zn  are  calculated  and  compared  with 
experiment.  The  temperature  dependence  of  these  coefficients  is  also 
analyzed  in  terms  of  changing  Fermi  surface-Brillouin  Zone  wall 
interactions.  In  PART  IV,  the  anisotropy  in  thermoelectric  power 
coefficients  is  calculated.  Finally,  a discussion  of  all  the  results 
presented  in  PARIS  II-IV  is  given  in  PART  V. 
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PART  II 


DRIVING  FORCE  FOR  ELECTROMIGRATION 

A.  General  Theory 

The  driving  force  on  a diffusing  ion  in  the  presence  of  high 
electric  current  derives  from  two  possible  sources.  The  first  is  due 
to  the  electrostatic  interaction  of  the  ion  with  the  external  field  and 
the  second  is  due  to  the  interaction  with  the  charge  carriers.  In  this 
treatment,  we  split  the  total  effective  charge,  Z*,  into 

z*  - Ze  + Zw  , (2.1) 

where  Z#  is  the  effective  electrostatic  charge  and  Zy  is  the  momentum 
transfer  or  electron  wind  charge.  (It  should  be  pointed  out  that 
although  this  separation  has  been  the  starting  point  of  most  discussions 
of  eleetroniigration,  it  is  not  essential.  Some  recent  investigators^”^1* 
have  preferred  to  deal  directly  with  the  total  driving  force  from  first 
principles  via  a linear  response  formalism.)  In  the  next  section,  we 
present  a few  remarks  concerning  Zg,  and  then  Zy  will  be  considered  in 
detail. 

1.  Electrostatic  Force 

A-  controversial  question  in  the  theozy  concerns  the  extent 
to  which  the  direct  force  of  the  applied  field  is  screened  by  the 
conduction  electron  gas.  In  the  models  of  Huntington^0  and  others^1***1 
there  is  no  appreciable  screening  of  the  external  field  and  the  full 
electrostatic  force  ZeE,  where  Z is  the  ionic  valence,  is  felt.  On 
the  other  hand,  Bosvleux  and  Prledel  7 maintain  that  for  interstitial 
ions  the  direct  field  force  is  totally  screened  out.  Recently,  however, 
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Shan  claims  to  have  shown  that  the  Bosvieux-Friedel  argument  Is  wrong, 

basing  his  results  on  a quantum  field  theoretic  treatment  of  the 
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problem.  Landauer  argues  that  significant  local  screening  may  be 
present,  tracing  the  source  to  what  he  calls  "spatial  conductivity 
modulation"^  and  "residual  resistivity  dipoles.  ^ i8 

3*4.  66 

generally  conceded  that  the  effects  suggested  by  Landauer  are  real,  * 
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the  magnitude  of  these  effects  remains  uncertain.  Bibby  and  Youdelis 

attempted  to  answer  the  question  experimentally,  looking  at  the 

electrcralgration  of  carbon  and  nitrogen  in  the  Hall  field  perpendicular 

to  the  current  direction.  These  authors  concluded  that  the  full 

electrostatic  force  corresponding  to  the  nominal  ionic  valence  was 
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present.  A recent,  as  yet  unpublished,  communication0^  suggests  that 
the  controversy  has  not  yet  been  fully  resolved.  Nevertheless,  in 
comparing  our  calculations  to  experiment,  we  will  subscribe  to  the  view 
that  the  full  electrostatic  force  is  exerted  by  the  external  field. 

2.  Force  on  a Diffusing  Atom  in  a Drifting  Gas  of  Bloch  Electrons 
We  calculate  here  the  effects  of  the  electric  current  on  an 
atom  diffusing  in  a crystal  lattice.  Let  ^ represent  the  coordinate  of 
the  diffusing  particle.  From  the  relationship  between  a coordinate  and 
its  canonically  conjugate  momentum  we  write 

Jj<P>  , (2.2) 

where  $ represents  the  exact  wave function  of  the  system  and  H is  the 

total  Hamiltonian.  We  proceed  now  in  the  spirit  of  the  Born-Oppenheimer 
6q 

approximation,  i.e. , the  migrating  ion  is  considered  to  be  a very  heavy 
classical  particle  Interacting  with  a quantum  mechanical  electron  gas. 

At  each  point  along  its  diffusion  path,  the  ion  is  taken  to  be  fixed  in 
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position  as  It  interacts  with  the  electrons.  For  such  a classical 
particle,  we  may  identify  ^ — as  a force  on  the  particle  and  write 
Eq.  (2.2)  as 


We  are  Interested  only  in  the  current  dependent  electron  driving  force 
on  the  ion,  so  in  Eq.  (2.3)  we  replace  H by  U,  the  electron  defect 
interaction  potential,  (in  the  case  of  interstitial  diffusion,  U 
may  be  a single  point  ion  potential.  For  diffusion  by  the  vacancy 
mechanism,  U must  Include  the  potential  of  the  jumping  ion  as  well  as 

OQ 

the  "semi-vacancies"  at  either  end  of  the  diffusion  path.)3  Distortion 
of  the  lattice  due  to  the  presence  of  the  defect  will  be  neglected. 

Our  starting  point  then  is 


t • -<f g|i>  • 

If  we  separate  jj  into  a product  of  electron  and  defect  parts  in  which 
the  position*'  of  the  ions  are  fixed,  and  take  the  electron  wavefunction 
to  be  a product  of  one  electron  functions,  we  arrive  at 

? - • (2. 

In  Eq.  (2.5),  f (k)  represents  the  electron  distribution  function  in  the 

presence  of  the  applied  electric  field.  It  of  course  must  be  perturbed 

0 

from  the  Ferml-Dirac  function  f in  order  to  reflect  the  presence  of 

current.  In  what  follows,  we  will  assume  that  f (k)  is  a known  function, 

having  been  determined  from  an  Independent  solution  of  the  Boltzmann 

equation  in  the  high  temperature  range.  In  our  numerical  computations, 

22 

the  actual  solution  of  Chan  and  Huntington  for  the  electron  mean  free 
patht  (see  APPENDIX  C)  will  be  used. 
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For  an  Isolated  interstitial,  take  U to  be  of  the  fora 
u(r-R),  r being  the  electron  coordinate,  *R  the  coordinate  of  the 
diffusing  ion,  keeping  in  mind  that  u denotes  an  unscreened  potential. 

(In  the  case  of  diffusion  by  the  vacancy  mechanism,  we  would  take 

U * u(r-S)  - u(r-R^)  - u(r-^)  , (2.6) 

R^  and  R^  denoting  the  positions  of  the  semi-vacancies.  We  will  return 
to  this  point  later.)  Fourier  expand  u as  follows: 


where 


^5 


u(q)  =J^J*u(x)e“^,xdjc  . 


(2.7a) 


(2.7b) 


C/I  is  the  crystal  volume.)  Using  Eq.  (2.7),  we  arrive  at 

. (2.6) 

Rote  that  Eq.  (2.6)  is  valid  in  either  the  vacancy  or  interstitial  case. 
Substituting  in  Eq.  (2.5)  gives 

F - ^ ^(^)iTu(q)e"iq*R^lfjjei^‘r|t^  . (2.9) 

We  now  have  to  calculate  the  one-electron  wavefunctions 
At  this  point  in  his  development,  Felt3  replaced  U everywhere  by  a 
screened  potential,  screened  as  it  would  be  in  the  absence  of  current 
in  a free  electron  gas.  He  then  calculated  the  electron  wavefunctions 
t using  perturbation  theory,  neglecting  any  other  contributions  to 
screening.  We  depart  from  that  approach  here  and  proceed  as  in 
Sorbello?0  treating  U as  a bare  unscreened  potential  in  Eq.  (2.5), 
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but  accounting  for  the  electron-electron  interactions  by  calculating 
the  t's  self -consistently  (in  the  Hartree  approximation)  to  first 
order  in  the  defect  potential.  Unlike  the  Sorbello  treatment,  ve  assume 
a weak,  local  potential  for  the  ions,  and  treat  the  case  of  Bloch 
electrons.  Our  zero-order  vavefunctions  therefore  will  have  a spatially 
varying  part  of  the  form 


*!<*>  - | a5®.1(5«)-7'  , 


where  the  sum  is,  in  general,  over  the  set  of  all  reciprocal  lattice 
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vectors.  Using  first  order  scattering  theory;  we  can  write  the 
wavefunction  tj*  arising  from  the  electron- defect  interaction  as 

l^£>  + . (2.11) 

In  Eq.  (2.11),  the  potential  W must  be  self- consistently  determined. 

The  Bloch  functions  are  the  self-consistent  eigenstates  of  the 
crystal  lattice  potential.  V therefore  includes  the  defect  potential, 
as  well  as  the  change  in  the  screening  potential  of  the  conduction 
electrons  due  to  the  presence  of  the  defect.  In  order  to  determine  tne 
latter,  we_need  an  expression  for  the  total  conduction  electron  charge 


density,  which  may  be  written 


i(?)  - If £)<£♦* 


(2.12) 


Fourier  transforming  Eq.  (2.12)  gives 


■>(?)  (2.13) 


Ve  wish  to  obtain  expressions  in  terms  of  the  zero-order  vavefunctions ; 
accordingly,  ve  use  the  relation 


<*4  ■ 


to  rewrite  Eq.  (2.13)  in  the  form 

How,  the  perturbation  expression  in  Eq.  (2.U)  can  be  used  to  evaluate 
the  matrix  elements  in  Eq.  (2.15).  The  results  are 


and 


<«*!*> -<W*>  t9i . t(e) 


(2.16) 


+ f- 


6(5)  - + i a 


(2.17) 


where  the  orthogonality  of  the  Bloch  functions  fJL.WU> 

k k 

has  been  used.  With  Eqs.  (2.l6)  and  (2.17),  Eq.  (2.15)  becomes 

n($)  * ^f(k)<  0£le"i^,rl 

<^kle’i^*r|^5?>X  0?lWl 


A&®- 


eft)  - 1(?)  -la 


<^fe 


■iq«r 


€(k)  +ia 


(2.18) 


The  first  term  on  the  right  side  of  Eq.  (2.18)  represents  the  zero-order 
charge  density;  i.e. , the  charge  density  in  the  absence  of  the  impurity, 
depending  only  on  the  Bloch  wavefunctlons  (f> . We  need  an  expression  for 
the  change  in  charge  density  due  to  the  presence  of  the  defect,  which 
cones  from  the  remaining  two  terms.  f?»i n ng  this  charge  density  An, 
we  get  from  Eq.  (2.18)  the  expression 


I 
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The  function  f(5)  may  be  split  into  two  parts  as  discussed  in  PART  I; 


C#  | 


f(£)  = f°(k)  + g (?)  (2.20) 

where  f°  is  the  Fermi  distribution  function  and  g(k)  is  the  change 
giving  rise  to  the  current.  Making  use  of  Eq.  (2.20)  and  the  following 


relations : 


At)  - f°(-£)  , 

€(k)  = e(-l)  , 
e(t)  - -g(-t)  , 


«(5)  -t(!?)  + lot  c(5)  - e(Z) 


- irt&|€(1c)  - t(?)J  , 


Eq.  (2.19)  simplifies  to 


“<*>  *><*!»  I 


2f°(5) 


[€(*)  - €(5) 


2i«g(k)5 


[€(k)-6(^]^‘ 


Next  we  explicitly  divide  the  potential  W into  two  parts. 


V « U + AV„  , 

B 


(2.21) 


(2.22) 


where  &Vg  means  the  change  in  the  screening  potential  due  to  the  defect. 

AV  can  be  further  broken  up  into  AV^  + AVC,  where  the  superscript  0 

denotes  that  part  of  AV  independent  of  the  current,  and  AVC  is  the 

8 


,-»T  m 
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current  dependent  part.  With  this  notation, 

W = U + AV°  + AVC  . 

8 S 


(2.23) 


Now  write 

W°  = U + AVg  (2.24) 

so  that  W°  represents  the  defect  potential  screened  as  it  would  be  in 
the  zero  current  case,  (it  will  be  instructive  to  write  some  of  our 
expressions  in  terras  of  this  potential,  which  we  may  assume  to  be  known 
from  the  solution  of  the  zero-order  problem. ) 

Before  proceeding,  we  can  perform  an  intermediate  check  by 
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comparing  Eq.  (2.21)  with  an  expression  derived  by  Stoddart,  et.  al. 
by  an  alternate  formalism.  They  considered  the  charge  distribution 
around  a weak  scattering  impurity  in  a stationary  gas  of  Bloch  electrons. 
From  Eq.  (2.20),  we  get  for  the  change  in  charge  density  in  the  zero 


current  case 


which  agrees  with  their  equation  (3.10). 

Poisson's  equation  (transformed)  can  now  be  used  to  relate 


(2.25) 


An  and  AV  . Specifically, 


AV  (q)  * An(£) 

8 q 

By  combining  Eqs.  (2.26)  and  (2.21),  we  arrive  at  the  following  expression 
for  AVg(?)  : 


(2.26) 
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- 2.  <^rle‘iq'r  I ^><^IW°  | ^g(^)B[c(t)-C(^J]  . (2.27) 

Q 4»  iC)K 

(in  the  second  tern  on  the  right  of  Eq.  (2.27),  W°  rather  than  W appears 

6ince  the  term  involving  AV°  and  g(k)  contains  the  product  of  two  small 

6 

quantities.)  From  Eq.  (2.27)  we  can  write 
AV°(q)  = W°(q)  - U(q) 

2 


Si 


y.  <*fki?  *1  h-><h-\ «°  I 


(2.28) 


Eq.  (2.28)  relates  the  unscreened  defect  potential  U to  the  defect 


potential  screened  in  the  zero  current  case  (W  ).  After  considerable 
algebra,  Eq.  (2.28)  can  be  generalized  as  follows: 


• 

• 

• 

uft-S) 

U(q) 

= 

U(q+3) 

• 

• 

I 

• • 

I 

w°(?-S) 

w°(® 

H°re4) 

• 

• 

• 

mm 

(2.29) 


'fc  here  can  be  interpreted  as  a "dielectric  matrix"  in  analogy  with  the 
free  electron  dielectric  function  €(q).  Its  rows  and  columns  are 
labeled  by  reciprocal  lattice  vectors  and  its  elements  are  given  by 


«(W,?4)  • sjj  -f  - bsL  X 2f°W  , 

* €(£)  - C(t+q+L) 


(2.30) 


we 


where  L,  G,  and  H are  reciprocal  lattice  vectors.  In  Eq.  (2.30), 
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have  used  the  notation  of  Sham  and  Ziman.  In  this  notation,  Eq.  (2.29) 
becomes 


€(q,q+G)W°(q+G) 


(2.296) 


For  notational  simplicity,  we  write  Eq.  (2.29)  in  matrix-vector  form: 

sc**  ^ 

U(q)  = 6 W°(q)  . (2.31) 

In  this  form,  Eq.  (2.31)  formally  resembles  the  scalar  equation 

U(q)  * £(q)w°(£)  which  relates  the  screened  and  unscreened  potentials 
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in  the  free  electron  case.  The  relationship  between  these  potentials 
in  the  presence  of  the  crystal  lattice  (Bloch  electrons)  is,  of  course, 
much  more  complicated. 

In  order  to  proceed,  we  assume  that  Eqs.  (2.29)  can  be 
inverted  and  solved  for  the  potentials  W^(^Hl),  and  formally  write 

W°0J)  * €*^(5)  . (2.32) 

To  further  simplify  the  later  discussion,  note  that  £ may  be  written 
as  follows: 

- SS  OC  aiS 

6 = I - M , (2.33) 

* » 

where  I is  the  identity  matrix,  and  D is  a diagonal  matrix  whose 
elements  are  given  by 


(2.34) 

The  elements  of  A follow  from  Eq.  (2.30).  We  now  need  an  expression 
for  AV^(q),  which  we  get  from  Eq.  (2.27): 


D(q+G,q+G) 


4*e£ 


- 


We4 


-iq«r 


e(£)  - e(£) 


(2.35) 


- ^ X <^|e-I?'7|  %><*y  W°|  )i?>8(')6 C«(E)-e(ta  . 

q Ji.  k,k 

Writing  Eq.  (2.35)  in  matrix  fora  in  analogy  with  Eq.  (2.29),  and.  using 
Eq.  (2.33),  gives: 

AV®(q)  = ~A  AV°(q)  - 2*i  DB  I?(q)  (2.36) 

or 

AV®(q)  = -2«i6  D B W°(q)  , (2.37) 

where  the  matrix  B is  given  by 

* X s(I)6g<k).«(S.M3<^|.-1'j45>-?|^Et5tr> 

• <^tl'1(54)'tl  h>  ■ <2-38> 

Eq.  (2.37)  is  the  central  result  of  this  section;  it  represents  the 
distortion  of  the  screening  cloud  around  the  defect  in  the  pi'eseuce 
of  the  current.  It  is  straightforward  tc  show  that  in  the  free  electron 
limit,  Eq.  (2.37)  goes  over  into 

AV®(q)  l2«i^W0(S)  ^G(k)&[€(k)-C(k+q)]  , (2.39) 
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which  agrees  with  the  expression  derived  by  Sorbello. 

We  are  now  ready  to  return  to  the  expression  for  the  force, 

Eq.  (2.9),  substituting  the  results  of  the  screening  calculation.  First, 
rewrite  Eq.  (2.9)  replacing  "q  by  in  the  summation: 

F - ^-iq  u^e1"**^;  f . (2.40) 

The  sxra  over  k in  Eq.  (2.4o)  is  Just  Jln(q),  which  was  written  out  in 
Eq.  (2.l8).  The  first  term  is  Just  the  charge  density  in  the  absence 
of  the  defect  (arising  from  the  unperturbed  Bloch  wavef unctions ) and 
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Is  not  of  interest  here.  We  also  get  a contribution  from  An°(q) 
which  is  just  the  change  in  charge  density  in  the  zero  current  case, 
evaluated  in  Eq.  (2.25).  Thi6  part  of  the  force  i6  not  of  interest 
either;  we  want  the  current  dependent  part.  This  must  cane  from  the 
current  dependent  change  in  the  charge  density,  which  we  can  get  from 
Eq.  (2.37): 

AnC(q)  = hV°(q) 

4*e  8 

- - 2%^-  e'1(q,q»'S)D(?+K,a*t)B(q-.t,q*M)W0(q-.M) 

4*e  K,L,M 


“ €1(^»5+L)B(q+L,^+M)W°(q+M) 

W°(q+M)  may  be  expressed  in  terms  of  the  unscreened  potential  as 
follows : 

W°(q+M)  * ^ €_1(q+M,5-£)u($*£) 

Also  we  note  that 


(2.41) 


(2.42) 


€(q^q-t-L)  = £ (q+L,q) 

|q+L| 


(2.43) 


which  follows  from  the  definition  in  Eq.  (2.30).  Eq.  (2.43)  implies 

9 l 

that  the  same  property  holds  true  for  b ",  which  is  proved  as  follows : 
€(q+M,q+<$)  fc_1(q+S  T+L)  » 6-j  £ 

^ ^ g €*(q+$»q+ft)  €~1(q-*$,q+t)  « ft*  ? 

• %l  lWlae*1*(WW) 
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4 

I 


l 

I 

I 


Summing  both  sides  over  M gives 

iSl Vx(f.fcO  - li+tf 

Using  Eqs.  (2.42)  and  (2.44)  in  (2.4l)  yields 


(2.44) 


An  (q)  * -2*i  (q+L, T+m)  £,  1 (q+L,q)  £," 1 ( q+ M, q+ G ) U ( q+G ) (2.4?) 

Mow  Eq.  (2.45)  may  be  substituted  in  Eq.  (2.40)  to  give  the  complete 
expression  for  the  current  dependent  force: 

* " -2«  ^^e  q B(q+L,q+M)  ^(q+M,q+S)U(q+G)  € ^ (q+£,q)u*(q)  . 

* ' ,M  (2.46) 

Eq.  (2.46)  can  be  made  more  explicit  by  substituting  for*  from  its 
definition.  Also,  we  can  make  the  complete  R dependence  of  the 
expression  explicit  by  substituting  for  U(q+£),  the  defect  potential. 

In  our  later  computations,  we  will  be  interested  in  self-diffusion  in 
Zn  single  crystals  (via  the  vacancy  mechanism).  Therefore,  for  U(q-Kj) 
we  use 

»($«)  - - e-‘(q«)-a>]  , (2.i,7) 

which  follows  from  Eqs.  (2.6)  and  (2.7).  Making  these  two  substitutions 
in  Eq.  (2,46)  results  in 

•u*(?)S(q,0,R)  , (2.48) 

where  we  have  collected  the  R dependence  of  the  expression  into  the 
ten 

S(q,0,R)  - e1^-1^)-*-  - *-*(**)•£,]  . (2.49) 


l 
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For  purposes  of  comparison,  it  is  instructive  to  consider 

the  form  of  Eq.  (2.48)  in  the  free  electron  limit.  Taking  the  matrix 

a 

€ over  into  the  scalar  function  • , and  replacing  the  Bloch  functions 
by  plane  waves  we  get 

F - - 1*2  g(k)fr(k'-k)|w°(?'-k)|  2S(k'-k,0,R)6[€(5)-€(^|].  (2.50) 

(FE)  h £,k'  *■  J 

We  take  g(k)  from  Eq.  (1.2),  T having  been  determined  from  solving  the 

Boltzmann  equation.  Eq.  (2.50)  is  apparently  the  sum  of  momentum 

transfers  f»(k-k')  for  transitions  between  states  on  the  Fermi  surface, 

multiplied  by  the  corresponding  "Golden  Rule"  transition  probability, 

calculated  using  a screened  potential  w°  = uyj|  . For  an  isolated  point 

defect,  it  would  be  exactly  that.  For  the  vacancy  case,  the  formula 

is  complicated  by  the  appearance  of  the  factor  S,  which  accounts  for 

the  more  complex  defect  structure  involved  in  vacancy  diffusion. 

Eq.  (2.48),  however,  is  not  subject  to  such  a simple 
interpretation.  Self-consistent  screening  of  the  potentials  in  the 
Bloch  electron  gas  is  accounted  for  formally  by  the  appearance  of 
elements  of  the  inverse  dielectric  matrix.  In  simple  terns,  the 
screening  cloud  surrounding  the  diffusing  particle  is  distorted  by  the 
presence  of  the  lattice,  even  the  zero  current  case,  its  distortion 
depending  on  the  position  R of  the  particle  relative  to  the  perfect 
lattice.  This  last  point  is  worth  some  discussion,  since  in  expressions 
such  as  Eq.  (2.29a),  the  R dependence  of  the  screening  is  not  explicit, 
and  perhaps  not  clear.  (Assume  for  the  moment  an  Isolated  point  defect 
to  simplify  the  writing.)  Recalling  Eq.  (2.7),  Eq.  (2.29a)  becomes 

u<J>e-i?-*.2'e(t^)w0<*S)e-1(*S)-»  , (2.51) 
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and  the  R dependence  of  the  screening  Is  explicit.  Another  point  of 
possible  confusion  remains.  The  phase  factors  In  Eq.  (2. 51)  appear  to 
depend  on  the  choice  of  origin,  which  should  otherwise  be  arbitrary. 
However,  from  Eq.  (2.30),  the  elements  of  C also  depend  on  choice  of 
origin  Implicitly  through  the  Bloch  wavefunctions  . In  fact,  starting 
from  Eq.  (2.10),  It  Is  easily  shown  that  a change  of  origin  introduces 
an  exactly  compensating  phase  factor  In  the  elements  of  7L  . Once  the 
origin  of  coordinates  Is  fixed,  this  arbitrariness  In  phase  is  removed 

and  C”  , and  simultaneously  from  Eq.  (2. 51) . 

The  factor  in  Eq.  (2.1*8)  also  deserves  some  comment.  It 
apparently  plays  a role  analogous  to  t>(£-  ?)  , or  momentum  transfer,  in 
the  free  electron  case  (Eq.  (2.50)).  From  the  form  of  the  matrix 
elements  in  Eq.  (2.1*8),  we  see  that  In  the  summation  q is  restricted  to 
take  on  values  of  the  form 

^>»T?/-k+*N  , (2.52) 

■# 

where  N is  a reciprocal  lattice  vector.  In  terms  of  plane  waves,  this 
is  suggestive  of  simultaneous  scattering  by  the  defect  and  Bragg 
reflection-by  the  crystal  lattice,  resulting  in  a momentum  transfer 
(k-k'  -N).  Strictly  speaking,  however,  we  know  that  the  electron  wind 
force  cannot  be  analyzed  exactly  in  terms  of  momentum  transfer  on  a 
per  collision  basis,  since  the  Bloch  wavefunctions  are  not  momentum 
eigenstates. 

B.  Application  to  Zinc 


1.  The  NFE  Model 

In  applying  Eq.  (2.48),  by  far  the  most  serious  complications 
arise  in  attempting  to  calculate  ? and  * 


Each  element  of  €.  involves 
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I 

I 


I 

I 


I 


I 


stns  over  the  Bloch  states  and  energies  as  given  by  Eq.  (2.30),  and 
? has  been  the  topic  of  considerable  general  discussion  in  the 
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literature.  Little  attention  has  been  paid  to  actual  numerical 
computations,  particularly  of  off-diagonal  terns.  Some  approximate 


calculations  relevant  to  insulators  and  semiconductors  are  discussed 
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by  Slnha  in  his  review  article.  For  metals,  directional  effects  in 


screening  have  been  largely  Ignored.  In  applying  our  NFE  model  to  Zn, 


we  will  follow  the  latter  course  and  replace  € with  the  free  electron 
scalar  form.  We  will,  however,  account  for  the  Bloch  nature  of  the 
xero-order  states  explicit  in  the  matrix  elements  in  Eq.  (2.43)  with 
a 2-OIV  model.  We  should  warn  that  the  consistency  of  this  procedure 
Is  not  clear;  the  directional  effects  being  ignored  are  not  estimable 
in  any  onvlous  way.  Arguments  that  they  are  small  for  a NFE  metal 
would  seem  to  apply  equally  well  to  the  band  structure  effects  being 
retained.  We  will  see  that  those  which  are  included  in  the  computations 
do  affect  the  quantitative  results  without  significantly  altering  the 
overall  qualitative  picture. 

in  Eq.  (2.43),  the  above  approximation  results  in  the  expression 


• . (2.53) 


w°  being  a screened  potential  (w°  » u/t  ).  This  formula  is  similar  to 


those  derived  by  Feit^®  and  Fiks^  except  that  It  contains  "cross  terms" 


of  the  form  w®(^V?5)w0#($),  Such  terms  were  absent  from  the  above 
mentioned  formulas  because  the  authors  both  assumed  a plane  wave  state 
for  the  diffusing  particle,  resulting  In  a sort  of  "pseudomomentum 
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conservation"  condition.  As  pointed  out  earlier,  the  diffusing  atom 
is  treated  here  as  a heavy,  classical  particle.  The  inclusion  of 
quantum  effects  in  the  diffusion  process  is  beyond  the  scope  of  this 
treatment. 


In  APPENDIX  A,  we  describe  the  semi -empirical  form  factor 
F(q)  which  will  be  substituted  for  the  spherically  symmetric  ionic 
potential  w®  in  Eq.  (2.53).  F(q)  will  also  be  used  to  construct  our 
Fermi  surface  model  (see  APPENDIX  B),  approximated  by  considering 
distortions  arising  from  one  pair  of  Brlllouln  zone  walls  at  a time. 
The  two  plane  wave  wavefunctions  consistent  with  this  model  have  the 
standard  form 


a(n)eik*r  ♦ 


^(n)e 


i(k+K)*r 


(2.5M 


where  a and  P are  given  by 

“(“>  ■ s[1  + 1’’-“'/)  (r-„)2  * 7 lA 
. t J|i . ,r.„/y  (r^  Jl/2 


(2.55s) 

(2.55b) 


and  g is  the  parameter  which  determines  the  band  gap  splitting  (see 
APPENDIX  B).  In  Eqs.  (2.55),  7*K/2ky  and  u-k^/kj,  as  in  Fig.  B.2.  The 
sign  of  p is  the  same  as  F(K)  for  n>7  and  reverses  for  n<7.  This 
form  for  the  wavefunctions  allows  explicit  evaluation  of  the  matrix 
elements  in  Eq.  (2.53). 

In  APPENDIX  C,  the  Chan,  Huntington  solution  of  the  Boltzmann 
equation  is  briefly  discussed  and  their  findings  are  summarized.  We  take 


«(?)  - ef.t(k)6[K?)-^3  , (2.56) 


where 


X.,c  ■ r«,0(k)  ' 


(2.57) 


the  subscripts  referring  to  the  direction  of  the  applied  electric  field, 

either  in  the  basal  plane  (a)  or  along  the  hexagonal  axis  (c).  The 

k dependence  of  T'  in  Eq.  (2.57)  turns  out  to  be  weak  (see  Fig.  C.2); 

therefore,  we  replace  T*  (k)  by  the  appropriate  average  also  displayed 

a,c 

In  Fig.  C.2.  This  was  the  procedure  followed  by  Chan  and  Huntington 
in  their  conductivity  calculations.  The  relevant  numerical  values  are 

<T>  = A7  x lO"1**  sec.  and  <r>  = .28  x 10-1^  sec.  The  k 

a c 

dependence  of  the  velocity  factor  in  Eq.  (2.57)  is  Included  as  prescribed 
by  our  2-OIV  model  of  the  Fermi  surface.  Fran  APPENDIX  B,  we  have 

V1*)  - -it-  C1  - 


»a(m) 


where  the  subscripts  relate  to  the  reciprocal  lattice  vector  which 
defines  the  axis  of  cylindrical  symmetry  for  each  particular  region. 
The  a and  c components  of  ~v  are  then  obtained  from  v#  and  ^ by  a 
straightforward  coordinate  transformation. 

We  turn  our  attention  now  to  the dependent  factor  in  Eq. 
(2.53).  Denote  the  direction  of  the  applied  field  by  the  subscript  i 
(i  being  either  a or  c in  our  calculations).  The  effective  force  for 
atomic  migration;  l.e. , Z^eE^  , is  obtained  from  the  R dependent  force 
by  averaging  over  the  diffusion  path.  That  is, 

rn 

Z^eE1  - -i-  ^F-d$  , (2. 

where  jf  • Rg  - is  the  vector  which  defines  the  diffusion  Jump,  and 
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is  its  canpoaent  la  the  1 direction.  This  result  follows  simply 
from  relating  the  increase  (decrease)  in  the  flow  of  atoms  in  a 
particular  direction  to  the  lowering  (raising)  of  the  potential  barrier 
for  diffusion.  (For  electromigration  in  the  a direction  in  hep  Zn, 

Eq.  (2.58)  is  not  precisely  correct.  The  relation  is  complicated  by 
the  fact  that  two  types  of  crystal 1 ographically  different  diffusion 
Juaps,  one  entirely  in  the  basal  plane,  and  one  plane  to  plane  (primarily 
In  the  c direction),  contribute  to  the  diffusion  coefficient  and  to  the 
effective  charge.  We  will  ignore  this  latter  contribution  when 
calculating  Z*.  The  estimated  error  involved  in  this  approximation 
Is  less  than  5 per  cent.)  Take  the  origin  of  coordinates  at  the  point 
R * R^«  Using  Eqs.  (2.53),  (2.56),  and  (2.57)  in  (2.58)  gives 


Zj  - -2*  C«(?)-«p]5  [«(^)-€(^l^  (^ 


(2.59) 


JL  cos  d 


i,i 


where  we  have  Introduced  the  notation  JL^*  J?cos  0^  Of  course,  if 
the  applied  field  is  in  the  direction  of  the  diffusion  Jump,  then 
cos  0^  1.  We  need  only  consider  the  real  part  of 

the  integral  in  Eq.  (2.59)  since  all  of  the  other  quantities  are  real 
in  our  model.  This  can  be  explicitly  evaluated  using  the  following 
result: 


*#> 

A'T) 


LtJf  ! 


at 


(2.60) 
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Eq.  (2.60),  together  with  the  definition  of  S((f,M,R^  in  Eq.  (2.49), 
results  in 


Re 


J^q,M,R)dRj  = S(q,M ,?) 

_ sln(M.Jf)  _ sin(*.J?)  - sln(i?*jf)  + 


^ — = sijJfrfojJ  > (2>6l) 

M*J  q-I 

-»  -*  , , . 

In  the  free  electron  limit,  M =*  0 and  Eq.  (2.6l)  reduces  to 

t&3jt)  ■ i - ?zl201 

q*X 

as  it  should?1* 

2.  Calculations  and  Results 

For  the  purpose  of  numerical  computation,  the  sums  over  k# 
and  k in  Eq.  (2.59)  can  he  transformed  into  integrals  in  the  usual 
way;  the  delta  functions  resulting  in  two  surface  Integrals  over  the 
Fermi  surface.  The  result  is 

F.*S.  F.S.  Wcos^ 


•F(q)F(|q*Ml)s(q,M,jO— ~ ~r 

v(5)  v(J) 


(2.62) 


where  C ^/l6*5f|2.  Cfl0  is  the  atomic  volune. ) The  explicit  fora 
of  the  matrix  elements  can  he  seen  by  substituting  from  Eq.  (2.10)  for 
the  Bloch  functions  The  result  is 

(r-^'4) 

a3^k^a?'  (&•%('*) 


■f  Jr ivi(!) 


*i-c 

F«S«  F«S« 


cos 


0, 


. (2.63) 

v(?)  t(S5 


- *.i*r ' 
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Using  the  2-OPW  vavefunctions  defined  by  Eqs.  (2.54)  and  (2.55),  the 
four-fold  sun  in  Eq.  (2.63)  over  reciprocal  lattice  vectors  gives 
rise  to  16  possible  terms  involving  appropriate  products  of  the  functions 
°(m),  0(h),  g(u#),  and  0(u#).  Each  of  the  sixteen  integrations  must 
be  performed  separately  in  the  full  calculation. 

In  performing  the  actual  integrations  using  our  model  Fermi 
surface,  cylindrical  symmetry  in  each  region  around  the  associated 
reciprocal  lattice  vector  (Fig.  B.2)  allows  us  to  write 

“ -7kxd*dV|)  “j^dk,,  -^d^du  . (2.64) 

Integration  over  each  region  is  then  carried  out  over  the  azimuthal 
angle  from  0 to  2«,  and  over  the  appropriate  limits  of  the  variable 
M*  The  Integrations  were  done  numerically  using  a four-dimensional 
trapezoidal  rule.  The  interval  sizes  chosen  for  the  numerical  integration 
correspond  to  dividing  the  m f)  and  \x*  regions  each  into  1728  equal 
area  subregions.  Attempts  to  refine  the  mesh  further  were  found  to  be 
impractical  in  terms  of  computer  time.  The  integration  procedure  was 
checked  using  a similar  triad,  function  and  the  known  result  was  reproduced 
to  within  a per  cent.  Rigorous  bounds  on  the  error  for  the  actual 
integrand  would  be  difficult  to  estimate.  We  feel  confident,  however, 
that  the  error  was  no  more  than  a few  per  cent,  certainly  insignificant 
considering  the  approximate  nature  of  the  overall  calculation. 

The  results  of  the  calculation  are  displayed  in  Table  2.1. 

The  calculated  values  represent  the  electron  wind  contribution  to 
the  effective  charge;  in  order  to  compare  with  the  experimentally 
determined  Z#,  the  electrostatic  contribution  (in  our  model,  Zft*  +2) 
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must  be  added.  The  calculations  were  carried  out  in  stages , each 

feature  of  the  model  being  introduced  separately  to  assess  its 

quantitative  effect.  The  first  calculation  was  performed  using  a 

spherical  Fermi  surface  and  1-CIV  wavefunctions  for  the  electrons. 

In  subsequent  calculations,  distortion  of  the  Fermi  surface  and  2-OIV 

39 

wavefunctions  were  added.  The  experimental  results  of  Routbort  ' 
are  also  displayed  in  Table  2.1.  a complete  discussion  of  the  results 
is  presented  in  PART  V. 


PART  III 


lcw-field  hall  effect 


A.  Hall  Effect  in  Anisotropic  Materials 
1.  The  Boltzmann  Equation 

The  Boltzmann  equation  in  the  presence  of  tooth  an  electric  field  "e 

■# 

and  a magnetic  field  H can  toe  written 

v(k)  - ^ [$(?)  x H • =^[jl(k)  -.X.(k}J  W(k,k*)  d3k' , (3.1) 

where  ^(k)  is  defined  toy 


f(k)  - f°(k)  = -e  He  ._X.(k)l 

L -*  b€ 


(3.2) 


We  are  assuming  here  quasi-elastic  scattering  between  electrons  on  the 
Fermi  surface;  the  validity  of  this  assumption  traces  to  the  fact  that 
the  energy  carried  by  the  phonon  responsible  for  electronic  transitions 
in  the  high  temperature  regime  is  very  small  compared  with  the  Fermi 
energy  carried  by  the  electron.  W(1t,T)  in  Eq.  (3.1j  represents  the 
scattering  probability  when  the  state  x is  known  to  toe  occupied  and  kx 
to  be  unoccupied.  Accordingly,  W possesses  the  property 

W(k,k}  = wfok)  . ^ (3.3) 

Price^  suggested  expanding  f()s)  (or  equivalently,  A(k)  ) 
in  powers  of  the  magnetic  field;  i.e. , 

.A.  * A + M + *N  + . . . , (3.*0 

is  independent  of  H,  M proportional  to  H,  etc.  Putting  Eq. 

(3.^)  in  Eq.  (3*1),  and  equating  terms  with  like  powers  of  H gives 

v « - T)  W d3k  , (3.5*) 
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-si 

1*- 

* 

j (M  - V)  W dV  , 

(3.5b) 

e I 

’S' 

[t.h.5] 

1 

| M = 

% 

pN  - V)  W d3k  , 

(3.5c) 

etc. 

In  calculating  the  Hall  effect  to  terms  linear  in  the  magnetic  field, 
only  Eqs.  (3.5a)  and  (3.5b)  need  be  considered.  It  appears  then  that 
calculating  the  linear  Hall  effect  involves  solving  two  successive 

integral  equations;  i.e.,  first  solving  Eq.  (3.5a)  fort  as  in  the 
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conductivity  analysis  of  Chan  and  Huntington,  then  substituting  this 
result  in  Eq.  (3.5b)  and  solving  for  M.  However,  Price  has  shown  that 
in  reality  only  Eq.  (3.5a)  need  be  solved  (for  the  linear  Hall  effect). 
The  reason  for  this  is  that  we  really  do  not  need  to  explicitly  evaluate 
the  vector  function  M.  We  only  need  to  calculate  its  effect  on  the 
current  density  j;  i. e. , we  need  to  calculate 


e2  £ 


■Ji®  5*  * 


d3k 


ft  It  ^3 


K v 


d€ 


d-'k 


(3.6) 


As  before,  we  assume  the  solution  of  Eq.  (3. 5a)  is  known  from  previous 

work.  This  allows  us  to  calculate  the  first  of  the  integrals  in  Eq. 

(3.6).  Staurting  with  Eq.  (3.5b),  we  arrive  at  an  expression  for  the 

second  integral  as  follows.  Multiplying  both  sides  of  Eq.  (3.5b)  by 
«♦  df®  -* 

A , integrating  over  k,  and  rearranging  on  the  left  yields 

dC 


r (♦  >»  Sf ® •* 

-f fronts 


k}  d3k  d3k 


d3k  d3k' 


(3.7) 


10*  { 


3^ 


Interchanging  k and  V in  the  second  integral,  and  using  Eq.  (3*3), 
we  get 

i f . n?  * ?t)  t $£  a . far  f(*  ^ . t » dVid’k 

tl  J 8«  J LJ  8C  b*  J 

- r*M  v — d3k  . (3 

J 84 

df°  df° 

(The  elastic  scattering  assumption  allowed  us  to  replace  — - by  — . ) 

8t'  8« 

The  right  hand  side  of  Eq.  (3.8)  is  seen  to  be  the  integral  required 
in  Eq.  (3.6).  Thus  the  current  for  the  linear  Hall  effect  is  given  by 

— E*  f(vx  7\)  X d\  t and  depends  only  on  the  properties  of  T 

tl  J d€ 


(3.8) 


— H • ( (vx  7x)  X — — dJk  , and  depends  only  on  the  properties  of  X. 

tl  J d€ 

calculated  from  Eq.  (3.5a). 

2.  The  Onsager  Reciprocal  Relations 

In  this  section  we  write  down  some  general  results  deriving 

from  the  Onsager  relations.  They  will  be  used  in  a later  section  in 
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applying  our  model  of  the  Fermi  surface.  The  Onsager  relation  for 
electrical  conduction  in  the  presence  of  a magnetic  field  states  that 

5 (H)~-  5(-H)  , (3.9) 

where  3 is- the  conductivity  tensor  and  ***  denotes  the  transpose.  In 
order  that  our  calculated  results  possess  the  required  reciprocity, 

Eqs.  (3.5)  and  (3*9)  require  that  fT  v d\  be  a symmetric  tensor, 

P n _o  ~ 8 4 

and  that  | M v — d^k  be  antisymmetric.  The  first  condition  follows 

J 8« 

immediately  from  Eqs.  (3.5a)  and  (3.3).  Using  Eq.  (3.S),  we  see  that 

the  second  condition  requires  that  H * J*  x ft)  t A be 

antisymmetric,  or  equivalently  that  K • /( dS  x 7 k)  T be  antisymmetric. 

* FsS» 

(Here  dS  is  the  vector  element  of  surface  area  on  the  Ferai  surface.) 

To  see  what  is  going  on  here,  we  write  as  a sum  of  dyads: 


(3.10) 


Using  Eq.  (3*10),  the  antisymmetry  condition  becomes 

• (dS  x bi  \)  = - jf\  » (dS  x ^ a.  b^) 

F*S.  *TsS* 

°r  r*  r + 1 

H •/  dS  x £(7x)  X + I J = 0 ( 

r «S#  1 

Bie  quantity  in  is  just  *$(&)  ; thus,  the  reciprocity  condition 


will  be  satisfied  if 


fi 


7(wT) 


If  the  Fermi  surface  is  a smooth,  continuous  closed  surface,  Eq.  (3.12) 
is  shown  to  be  true  by  dividing  the  surface  in  two  parts,  and  applying 
the  generalized  Stokes  Theorem  to  each  part.  The  two  line  integrals 
which  result  cancel  exactly.  For  a surface  cut  by  Brlllouin  zone  walls, 
Stokes  Theorem  gives 


fis  X 7(Sx)  = Zi* 

«/  r r> 


where  each-term  in  the  sum  is  a line  integral  over  the  boundary  curve 
of  an  open  portion  of  the  Fermi  surface.  The  sum  in  Eq.  (3.13)  is  seen 


to  be  zero  by  considering  contributions  from  curves  on  opposite  faces 
of  the  Brlllouin  zone;  since  they  are  displaced  by  a reciprocal  lattice 
vector,  symmetry  requires T to  be  the  same  for  corresponding  points, 
whereas  the  direction  of  integration  is  opposite  (as  dictated  by  proper 
application  of  Stokes  Theorem).  Thus  our  solutions  will  have  the 
required  Onsager  symmetry,  provided  we  use  a continuous,  smooth  Fermi 
surface.  Interrupted  only  by  Brlllouin  zone  discontinuities.  Our 
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model  surface  (APPENDIX  B)  does  not  meet  that  requirement,  however. 

This  point  will  be  discussed  in  section  III.B.3. 

3.  Expressions  for  the  Hall  Coefficients 

The  HalT  constant  as  customarily  defined  is  given  by  the 

ratio  of  the  transverse  electric  field  to  the  product  of  the  current 

density  (in  the  direction  of  the  applied  electric  field)  and  the  magnetic 

field.  For  example,  if  the  experimental  geometry  is  such  that  the 

current  is  constrained  to  flow  in  the  "y"  direction,  and  the  magnetic 

field  applied  in  the  "z"  direction,  then  the  Hall  coefficient  is  given 

by  R * ~2lL.  . (x,y,  and  z form  a right-handed  coordinate  system  and 

yz 

the  sign  of  R is  such  that  free  electron  conduction  produces  a negative 
Hall  coefficient. ) The  Onsager  relations  allow  at  most  three  Independent 
Hall  coefficients,  depending  only  on  the  direction  of  the  applied 


magnetic  field.  For  H 

in  the  x direction,  one  obtains  using  Eqs. 

(3.6) 

and  (3.9) 
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where  C * U«3ti/e2  and  4^=  >J^’ivi“v  is  ProPortl°nal  to  the  conductivity 
in  the  absence  of  a magnetic  field.  (All  of  the  integrations  in  Eqs. 
(3»l^)-(3.l6)  are  over  the  Fermi  surface.)  The  expressions  in  the 
numerators  of  Eqs.  (3.1*0-(3.l6)  can  be  evaluated  in  terms  of  the 
vector  function  X.  and  its  derivatives  by  making  use  of  Eq.  (3.6). 

Making  the  appropriate  substitutions,  ve  arrive  at  our  final  expressions 
for  the  three  independent  Hall  constants: 


r(ha) 


e AJJ  6kk 


f(y  ^ _ T x 2£ 

J J 4k  * ik.  J V 


vhere  e is  taken  to  be  the  (negative)  electronic  charge,  and  i,J,k 
represent  a permutation  of  the  variables  x,y,  and  z.  (The  two 
expressions  for  R(H^)  generated  by  Interchanging  J and  k in  Eq.  (3.17) 
are  of  course  equal  as  a result  of  the  Onsager  relation.) 

B.  Application  to  Zinc 
1.  General  Formulas 


As  before,  we  take  and  substitute  into  Eq.  (3.17) 


to  get 


r(ha)  - 


— - v.  -Jew 


1 - fan- kvjY 

J*  dr  dr 

The  second  Integral  In  Eq.  (3«l3)  contains  (v  x Vt)^  , and  thus  involves 


r 
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non-uniformities  in  7*  from  point  to  point  on  the  Fermi  surface  (since 

v is  perpendicular  to  the  surface) . For  R(H  ),  only  the  directional 

z 

derivatives  of  T'  along  orbits  parallel  to  the  basal  plane  enter. 

However,  the  Chan,  Huntington  calculations  ’ demonstrated  that  this 

variation  is  slight  in  general.  Justifying  its  neglect  here.  For 

8(1^)  (or  R(Hy)  ) the  situation  is  not  quite  so  clear.  In  these  cases 

the  variation  of  T with  u,  as  shown  in  Fig.  C.2,  is  involved.  The 

deviations  of  T from  <T>  still  appear  small,  though  not  so  clearly 

negligible  as  in  the  first  case.  From  a practical  point  of  view,  it 

is  doubtful  whether  the  calculated  data  are  accurate  enough  to  Justify 

Inferring  values  of  the  derivatives  of  on  a point  by  point  basis. 

nevertheless,  using  the  data  of  Fig.  C.2  and  a spherical  Fermi  surface 

corrected  for  the  missing  area  (see  Fig.  4.2),  a rough  estimate  of  the 

second  Integral  in  Eq.  (3- 18)  was  calculated.  The  calculation  proved 

to  be  a fairly  sensitive  one,  owing  to  close  cancellation  between  regions 

of  positive  and  negative  — . The  calculated  estimate  was  roughly 

20  per  cent  of  the  contribution  from  the  first  integral,  and  of  the 

same  sign.-  In  what  follows,  the  contributions  from  Vf*  are  not  included. 

For  Zn,  hexagonal  symmetry  reduces  the  number  of  independent 

Hall  constants  to  two,  one  with  the  magnetic  field  in  the  c direction, 

and  the  other  with  the  magnetic  field  in  the  basal  plane.  Choosing 

our  z-axls  along  the  hexad  symmetry  axis,  we  need  to  consider  R(H  ) 

z 

and  either  R(HX)  or  R(Hy).  We  make  the  following  definitions: 

Ig  ■ R(HX)  and  R^«  R(Hy).  From  Eq.  (3*18),  we  arrive  at  the  following 
expressions  (with  - <'T^>vi)  : 


39 


1 bv  bv 

/ (y  _JL  . v _Jhy  ds 

W* J x\yy 

* (Vi*.  ("v2dS 

(3.19a) 

J x v J y v 

- (***y) 

(3.19>) 

f(r  t li  . „ fli),  M 

„ * ’ 
A ' (V  32  . P i£ 

J X V J Z V 

(3.20a) 

- (*«♦*) 

(3.20b) 

where  the  notation  (x  ♦*  y ) means  Interchange  the  variables  x and  y in 
the  preceding  expression.  Eqs.  (3.19)  and  (3.20)  are  the  starting 
points  for  our  numerical  computations. 

In  performing  the  required  calculations  for  our  Fermi  surface 
model,  we  encounter  a difficulty  suggested  in  section  III. A. 2.  The 
Onsager  relations  require  that  our  model  yield  the  same  results 
whether  the  a or  b versions  of  the  above  equations  are  used.  By 
symmetry  considerations  SLlone,  it  is  easily  seen  that  Eqs.  (3.19a) 
and  (3.1S*>}_will  give  the  same  result.  However,  expressions  (3.20a) 
and  (3.20b)  do  not,  in  general,  give  the  same  result.  The  problem 
arises  because  we  have  Introduced  artificial  boundaries  in  our  model 
surface;  the  free  surface  area  is  divided  up  into  cyllndrically 
symmetric  regions  about  the  appropriate  reciprocal  lattice  vectors. 
There  is  a certain  amount  of  arbitrariness  in  the  relative  size  of 
these  regions.  We  can,  at  least  in  our  model  for  Zn,  remove  this 
arbitrariness  by  selecting  our  boundaries  so  that  the  Onsager  relations 
are  satisfied.  The  required  condition  from  Eq.  (3.13)  is,  in  general, 
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2 (*A) 


c c 


(3.21) 


The  particular  component  of  the  triad  of  interest  here  is  the  J 
component;  the  condition  we  need  to  satisfy  in  order  that  Eqs.  (3.20a) 
and  (3. 20b)  give  the  same  result  is 

dkvv  = 0 . (3.22) 


U- 


y x z 


The  sum  over  the  curves  C representing  intersections  with  Brillouin 
zones  gives  zero  (as  was  previously  argued),  which  leaves  the  line 
integrals  over  the  artificial  boundaries  introduced.  The  limits  of 
integration  in  each  i *gion  were  varied  (actually,  symmetry  and  the 
requirement  of  constancy  of  free  Fermi  surface  area  reduced  the 
nunber  of  variable  paLrameters  to  one)  until  the  sum  in  Eq.  (3.22) 
was  forced  to  zero.  This  determined  the  limits  of  integration  used 
in  the  final  calculations. 

2.  Calculations  at  High  Temperature  (6$0  K) 

In  calculating  the  Hall  constants  the  following  procedure 
was  used.  _ln  each  cylindrical  ly  symmetric  region,  the  tensor 

x V v)$  ^ was  calculated  with  the  z-axls  chosen  along 
the  axis  of  symmetry.  This  made  fa  diagonal  and  also  simplified 
the  expressions  for  v.  These  are  given  explicitly  in  APPENDIX  B. 

Then  a straightforward  similarity  transformation  was  used  to  transform 
the  result  to  a coordinate  system  with  z-axls  in  the  c direction.  The 
calculated  tensor  H • ^ (v  x Vv)T  proved  to  be  antisymmetric, 
which  verified  the  procedure  described  in  the  previous  section  for 
selecting  the  limits  of  Integration.  The  integrals  in  the  denominators 


of  Eqs.  (3.19)  and  (3.20)  were  calculated  In  a similar  manner. 

The  high  temperature  Fermi  surface  parameters  relevant  to 
our  model  are  shown  In  Table  B.l.  The  calculated  results.  In  units 
of  mJ/coul. , are 

R|l  ■ +8*2  * 10 

- -1.3  » 10"U  . 

By  comparison,  the  free  electron  value  for  Zn  is  about  -4.75  x 10“^. 

The  values  are  shown  in  Fig.  3.1,  along  with  several  experimentally 
determined  values. 

3.  Temperature  Dependence 

Of  course,  a complete  temperature  analysis  would  require 
resolving  the  Boltzmann  equation  for  ^(?)  at  another  temperature  or 
set  of  temperatures.  Short  of  this,  if  we  choose  another  temperature 
and  again  make  the  approximation  that  at  this  new 

temperature,  T*  again  disappears  from  our  expressions  for  R„  and 
R^.  These  expressions  depend  only  on  the  size  and  shape  of  the  Fermi 
surface  whose  temperature  dependence  can  be  handled  approximately 
within  our~model.  400  K was  chosen  as  the  temperature  for  a second 
set  of  calculations,  and  the  relevant  lattice  and  Fermi  surface 
parameters  were  recalculated.  A comparison  between  these  values  and 
those  given  in  APPENDIX  B for  650  K is  given  in  Table  3.1.  The  values 
of  c and  a were  recalculated  using  the  thermal  expansion  data  of 

do 

OUder  and  Wallmark.  In  addition  to  changing  the  c/a  ratio,  increasing 
the  temperature  tends  to  decrease  the  Debye-Waller  factors®^  (e-2W)  and 
thus  decreases  the  band  gaps  responsible  for  Fermi  surface  distortion. 

The  results  of  the  calculations  at  400  K are  also  shown  in  Fig.  3.1. 
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The  available  experimental  data  for  the  Hall  effect  in  Zn 

single  crystals  extend  only  to  about  300  K.  However,  polycrystalline 
54  66 

data  * for  both  Zn  and  Cd  decrease  approximately  linearly  with 

temperature  in  the  region  beyond  300  K.  In  addition,  single  crystal 

87 

measurements  for  both  R)(  and  in  Cd'  ' show  a roughly  linear 
decrease  with  temperature  in  this  range;  all  of  which  suggests  a 
linear  extrapolation  of  the  experimental  data  in  Fig.  3*1  for  Zn.  Our 
calculated  values,  also  6hown  in  Fig.  3-l>  would  apparently  agree  quite 
well  wj t ay  reasonable  extrapolation. 


PART  IV 


THERMOELECTRIC  POWER 


A.  General  Theory 

In  the  presence  of  a temperature  gradient,  the  Boltzmann 
equation  for  the  conduction  electrons  becomes  (again  assuming  quasi- 
elastic  scattering) 


(-  *.[•?-  =J[e&  - e(*)]  dV  . (4.1) 

Extending  the  mean  free  path  formalism,  one  may  define  * by 


6(k)  * (-  ^)  [eE  - (€  -€p)fc]  . X . (4.2) 

Substitution  of  Eq.  (4.2)  in  (4.1)  results  in  the  simpler  form  of  the 
Boltzmann  equation. 


v(k)  = ^J[X(k)  - ‘£(ls)J  W(k,V)  . (4.3) 

In  anisotropic  materials,  the  absolute  thermoelectric  power  (TE1 ; is 
defined  by 

Ea  « Si(^T)i  , (4.4) 

where  the  subscript  i denotes  the  direction  of  the  temperature  gradient. 

In  terms  of  T,  the  resulting  expression  for  becomes 


where  6i(€)  is  given  by 

61(‘>  - j5i  • <‘-6> 
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In  Eq.  (4.6),  the  integration  is  over  a surface  of  constant  energy 
€.  Eq.  (4.5)  is  to  be  interpreted  as  follows.  It  means  that  one 
should  calculate  the  conductivity  as  a function  of  the  hypothetical 
height,  €.  , of  the  Fermi  level,  and  then  evaluate  the  derivative  of 
<S(€)  with  respect  to  6 at  the  actual  Fermi  energy,  € . The  variation 

r 

of  with  energy  will  depend  on  the  change  in  the  numbers  of  electrons 
available  for  conduction,  as  well  as  the  change  of  the  mean  free  path 
of  each  electron  as  its  energy  changes. 


With  TjV^  , Eq.  (4.6)  becomes 

2 
e 


JL  2 dS  - e2  . f 2 dS  . e2  „ 

where  = ^"v2  ~v  * 411(1  ^ ) denotes  the  appropriate  average  value. 
Substituting  Eq.  (4.7)  in  (4.5)  gives 


(M) 


, ^ r«p  \ + a<ri> " 

Ji  ' |^xi  +<rA>  ae 


(4.8) 


€=€. 


Eq.  (4.8)  is  the  form  in  which  the  TEP  is  often  analyzed.  Since  e 
is  negative,  the  factor  preceding  the  dimensionless  quantities  in 
brackets  is  negative  and  typically  of  the  order  of  one  microvolt  per 
degree  K.  The  first  term  in  the  brackets  in  Eq.  (4.6),  call  it  s^  , 
involves  only  the  geometry  of  the  Fermi  surface,  and  has  the  value 
3/2  in  the  FEM.  Ignoring  the  second  term  for  the  moment,  this  would 
predict  a negative  TEP.  In  his  analysis  of  the  noble  metals,  Ziman**1' 
showed  that  the  effect  of  Fermi  surface  anisotropy  was  to  reduce  s^  , 
but  not  enough  to  drive  it  negative  and  account  for  the  positive  values 
of  TEP  observed  experimentally.  The  experimental  results  were  only 
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understood  after  a careful  analysis  of  the  tem  involving  <?*>,  which 
we  shall  hereafter  call  Sg. 

B.  Application  to  Zinc 

1.  Fermi  Surface  Model 

In  our  TEP  calculations  for  Zn,  we  chose  a simplified  model 
Fermi  surface  illustrated  in  Fig.  4.1.  The  surface  is  spherical,  and 
1-OIW  electron  wavefunctions  are  used  throughout.  The  missing  area 
(APPENDIX  B)  was  accounted  for  approximately  by  multiplying  each  element 
of  surface  area  by  P(u),  where  P(u)  is  given  by 

w>1w 


and 


cos 


* 


m 


r101  • uc°eei0i 
sin  0, 


(4.10) 
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? ~L01 

being  the  angle  between  the  reciprocal  lattice  vectors  Q_  (see 
Eqs.  (4.15))  and  the  c-axis.  In  equations  (4.9)  and  (4.10),  the 
quantities  UT|nY,  and  will  of  course  depend  on  6 , the 

hypothetical  Fermi  level  in  Eq.  (4.8).  This  approximate  treatment  of 
the  missing  area  assumes  it  to  be  uniformly  distributed  in  the  ^ direction. 

2.  Solution  of  the  Boltzmann  Equation 

°2 

Following  Chan  and  Huntington7  we  write  * ^\v±  • from 
which  it  follows  that 

- v?) 


vt(k) 


j'fvjft)  - v^S)]  W(k,k}  d3k" 


(*».U) 


c-axis 


Fig.  4.1!  Fermi  Surface  Model  (TEP  Calculation) 
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Strictly  speaking,  Lq.  (4.11)  actually  represents  the  first  term  in 
an  iterative  expansion  for  T^.  However,  subsequent  terms  were  shewn 
to  be  very  small,  and  we  shall  therefore  use  Eq.  (4. 11)  for  in 

this  calculation.  What  is  needed  then  is  a solution  for  T and  T* 

a c 

as  a function  of  the  hypothetical  Fermi  level, € . in  the  actual 
calculation,  the  energy  derivatives  were  estimated  by  calculating 
ri<«“V)  and  T^(€=tp  + Afp),  with  A taken  to  be  .002. 

We  next  describe  an  explicit  fora  for  the  scattering  kernel 

W(k,$),  the  basic  transition  probability  from  state  k to  We  use 
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the  rigid  ion  model  for  the  electron  phonon  interaction,  along 
with  a local  pseudopotential  approximation.  The  optical  phonon 
branches  are  accounted  for  by  using  the  double  (or  Jones)  zone  (Fig. 
B.l)  for  the  phonon  cells.  Although  this  has  a complex  shape,  it  is 
not  too  far  from  spherical;  we  therefore  replace  the  actual  phonon 
cells  by  Debye  spheres.  The  experimental  phonon  dispersion  curves 
and  our  analytic  fits  to  these  curves  are  described  in  APPENDIX  C 
and  illustrated  in  Figs.  C.l  and  C.2  respectively. 

The  contribution  of  normal  phonon  processes  (i^)  to  the 
Integral  I (it)  is  given  by 


y?)  = C fas  (-q)i-V  F2(q)  , (4.12) 

J V 

■*/  -» 

where  q = k - k is  the  phonon  wave  vector,  F(q)  the  pseudopotential 
the 

V 


fora  factor,  <jd  t the  (longitudinal)  phonon  frequency,  and 
qi. 
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In  Eq.  (4.13),  kp  is  the  Fermi  vavevector,  n the  number  of  ions  per 
unit  cell  (2  for  Zn),  N the  number  of  cells  per  unit  volume,  and  M 
is  the  ionic  mass.  The  domain  of  Integration  over  k implicit  in  Eq. 

(4.12)  is  determined  by  the  intersection  of  the  Debye  and  Fermi  spheres, 
and  1b  illustrated  in  Fig.  4.2(a).  The  contribution  of  umklapp  processes 
Is  given  by 

yk)  = C J"ds'£(£ + '$)]1  J2(q,$)  F2(q+Q)  s|  , (4.l4) 

where  q + Q = k - k and  Q is  a reciprocal  lattice  vector.  is  the 
structure  factor  for  Q defined  in  the  usual  way  by 


n-1  -* 

e19-r„ 


where  r specifies  the  position  of  each  atom  within  the  unit  cell, 
m 


(^.15) 


J2(q,Q)  is  defined  by 


Ati)  = Bl:  . ilSjili 


(*•16) 
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e^  being  a unit  vector  in  the  q direction,  and  to  means  the  polarization 
is  transverse  and  lying  in  the  q-Q  plane.  The  geometry  for  urtklapp 
scattering  and  the  domain  of  integration  for  Eq.  (4.14)  are  schematically 
illustrated  in  Fig.  4.2(b). 

3.  Calculations 

In  computing  the  umklapp  integrals,  contributions  from  four 

types  of  RLV's  were  considered.  With  the  z direction  along  the  crystal 

c-axls,  the  components  of  the  RLV's  considered  are 
4* 


^ ^ (0,0,±1) 
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Debye  Cell 


(a)  Normal 


Fig.  4.2!  Geometry  for  Normal  and  Umklapp 
Scattering 
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01  these  32  HLV's,  18  give  non- zero  contributions.  The  surface 
Integrations  were  performed,  numerically,  and  the  area  elements  used 
were  roughly  1/8,000  of  the  total  Fermi  sphere  area. 

The  first  set  of  calculations  was  performed  using  the  model 
form  factor  described  in  APPENDIX  A (dotted  line  in  Fig.  A.l).  The 

calculated  values  of  T and  T*  are  plotted  vs.  4 in  Fig.  4.3.  (We 

a c 

should  point  out  that  a singularity  in  the  scattering  kernel  occurs 

near  4 = . 8 , as  1 reaches  the  Brillouin  zone.  This  type  of  singularity 

is  known  tc  be  a result  of  the  FEM;  it  does  not  occur  when  the  proper 

Bloch  wavefunctions  are  used.  The  plotted  results  represent  a smooth 

continuation  over  the  singularity. ) The  comparable  results  from  Ref. 

22  are  shown  in  Fig.  C.3.  Results  of  the  TEP  calculations  for  and 

S are  given  in  Table  4.1,  along  with  values  of  the  dimensionless 
c 

quantities'sj^  and  s^.  Results  of  a second  set  of  calculations  using 

the  semi-empirical  form  factor  (APPENDIX  A)  are  also  given  in  Table 

4.1.  Fig.  4.4  displays  the  theoretical  values  and  the  experimental 

results  of  Rowe  and  Schroeder.  ^ 

The  experimental  data  for  Zn  single  crystals  are  quite  sparse, 

extending  only  to  about  room  temperature.  As  in  the  Hall  effect,  a 

straight  line  extrapolation  to  high  temperature  is  again  suggested  by 
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the  available  polycrystalline  data.  Based  on  this  assunption,  our 
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Fig.  4.3!  Calculated  Values  of  t vs. /x. 
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Fig.  4.4'.  Thermopower  vs.  Temp. 


calculations  appear  to  agree  qualitatively  with  experiment.  The 
contributions  of  the  individual  terms  s^  and  will  be  discussed 
In  the  next  section. 
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PART  V 


DISCUSSION  AND  CONCLUSIONS 

This  work  has  dealt  with  the  calculation  of  the  high  temperature 

transport  parameters  of  single  crystal  Zn.  Zn  is  a divalent  metal,  but 

has  small  band  gaps;  consequently,  we  have  employed  a 2-OIV  NFEM 

(APPENDIX  B).  Fermi  surface  experiments  indicate  that  the  third  band 

"butterflies"  in  Zn  are  unoccupied.  This  is  consistent  with  our  model 

(missing  area)  which  predicts  that  they  remain  unoccupied,  even  at  high 

temperature.  We  have  also  used  a local  pseudopotential  approximation, 

which,  as  discussed  in  APPENDIX  A,  is  adequate  for  the  calculation  of 

transport  properties.  Such  a model  is  able  to  account  for  the  high 

22 

temperature  conductivity  of  Zn;  our  calculations  have  demonstrated 
its  applicability  to  a variety  of  other  transport  phenomena.  Cd  exhibits 
similar  transport  property  anisotropy,  so  that  much  of  the  discussion 
below  applies  qualitatively  to  Cd  as  well. 

A.  Driving  Force  for  Electromigration 

In  PART  II.A. , we  have  reconsidered  band  structure  effects  in 
electromigration  by  calculating  the  driving  force  on  a weak  scattering 
impurity  in  a drifting  gas  of  Bloch  electrons.  Our  approach  differs 
from  previous  treatments  in  that  we  explicitly  include  directional 
effects  in  conduction  electron  screening  via  a self-consistent  calculation 
of  the  scattered  electron  wavefunctlons.  As  a result,  our  formal 
expression  for  the  driving  force,  Eq.  (2.48),  depends  explicitly  on 

fa 

the  dielectric  matrix.  Although  we  have  introduced  € in  a formal  way, 
in  our  numerical  computations  we  have  not  gone  beyond  the  approximation 
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of  treating  £ as  a scalar.  In  this  limit,  our  formula  reduces  to  those 
of  Feit  and  Fiks,  except  for  the  restriction  they  impose  by  treating 
the  migrating  atom  as  if  it  were  in  a plane  wave  state.  In  the  limit 
of  vanishing  lattice  potential,  Eq.  (2.48)  agrees  with  the  result  of 
Sorbello. 

The  agreement  between  our  calculated  values  and  experiment 
must  be  considered  quite  good,  even  though  the  predicted  anisotropy 
is  somewhat  low.  There  is,  unfortunately,  a certain  amount  of  arbitrariness 
in  some  of  the  factors  which  contribute  to  the  computed  values  of  Z* 

First,  although  we  calculate  Z , in  order  to  compare  with  experiment 
we  assume  that  Zg  = +2,  corresponding  to  the  full  electrostatic  force 
being  present.  Although  we  feel  this  is  essentially  correct,  it  rerains 
a somewhat  controversial  question.  Second,  the  values  of  are 

taken  from  an  independent  solution  of  the  Boltzmann  equation.  However, 
since  they  are  approximately  consistent  with  the  experimental  conductivity, 
they  might  alternatively  be  considered  as  empirically  determined  parameters. 
In  fact,  closer  agreement  with  the  experimental  anisotropy  would  have 
been  obtained  if  the  relaxation  times  had  been  fitted  to  the  actual 
experimental  conductivity  values.  Finally,  the  calculated  values  of 
Zy  are  of  course  dependent  on  choice  of  pseudopotential  form  factor.  The 
use  of  the  model  form  factor  (dotted  line  in  Fig.  A.l),  for  example, 
yields  values  roughly  30-40  per  cent  greater  than  those  shown  in  Table  2.1. 

Despite  the  qualifications  discussed  in  the  previous  paragraph, 
several  unambiguous  conclusions  emerge  f ran  the  numerical  computations. 

Quite  clearly,  the  anisotropy  in  effective  charge  is  due  primarily  to 
the  anisotropy  in  relaxation  times.  Unlike  the  conductivity  calculation. 


vhere  Fermi  surface  distortion  compensates  the  anisotropic  relaxation 

times  quite  closely,  only  a slight  compensation  occurs  in  the  electron 

vind  force.  We  find  that  Z */Z,C  = 1.35,  compared  to  K'T  >/<T*>  = 1.65. 

vr  v a c 
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This  confirms  the  suggestion  of  Chan  and  Huntington  based  on  rough 
estimates  of  the  wind  forces. 

Table  2.1  also  points  out  the  magnitude  of  the  various  band 
structure  effects  which  we  have  considered.  The  effect  of  the  missing 
Fermi  surface  area  is  quite  drastic  since  the  wind  force  formula 
involves  two  integrations  over  the  Fermi  surface.  Another  result  of 
Fermi  surface  distortion  is  the  vanishing  normal  component  of  electron 
velocity  (and  mean  free  path)  at  a zone  boundary;  this  tends  to  reduce 
transport  perpendicular  to  the  zone  boundary.  There  is  also  the  fact 
that  Bloch  waves  are  not  momentum  eigenfunctions,  so  that  a given 
scattering  event  does  not  give  rise  to  a well  defined  momentum  transfer. 
In  terms  of  free  electrons,  an  electron  may  be  simultaneously  scattered 
by  the  defect  and  Bragg  reflected.  We  find  that  this  also  tends  to 
decrease  the  effective  charge. 

We  should  also  comment  on  the  effect  of  defect  structure  in 
the  case  of  vacancy  diffusion.  Although  the  full  calculations  were  not 
repeated  for  an  isolated  point  defect,  some  preliminary  calculations 
indicate  that  the  effective  charge  for  the  vacancy  case  is  roughly  30 
per  cent  lower  than  for  the  interstitial  case. 

B.  Hall  Effect 

At  650  K,  our  calculated  value  for  R ^ is  roughly  1.5  times 
the  free  electron  value  for  Zn  in  magnitude,  and  opposite  in  sign.  This 
large  positive  value  is  due  to  the  regions  of  negative  curvature  in  the 


"monster"  section  of  the  Fermi  surface,  i.e.,  to  the  region  between 
horizontal  (002)  Brillouin  zone  planes.  In  our  model,  this  region  is 
strongly  distorted  by  the  slanting  (101)  zone  planes.  Fig.  B.2 
indicates  the  reversal  of  curvature  below  the  BZ  intersections.  The 
free  electron-like  "lens"  (polar  caps),  on  the  other  hand,  has  positive 
curvature  and  contributes  toward  a negative  Hall  constant.  Compensation 
between  regions  of  opposite  curvature  result  in  a small,  negative  value 
for  R^ . 

As  discussed  in  PART  III,  our  numerical  values  depend  on  the 
somewhat  arbitrary  division  of  the  free  Fermi  surface  area  into  regions 
(Fig.  B.3).  Although  our  choice  of  boundaries  based  on  the  Cnsager 
symmetry  requirements  seems  physically  reasonable,  it  would  be  more 
satisfying  to  have  a model  without  this  degree  of  arbitrariness.  The 
Rj.  calculation  turned  out,  in  fact,  to  be  quite  sensitive  to  the  precise 
choice  of  boundaries  owing  to  close  cancellation  between  positive  and 
negative  contributions.  The  R()  calculation,  however,  proved  to  be 

nearly  independent  of  this  choice.  The  same  conclusion  can  be  drawn 

■* 

concerning- the  neglect  of  the  V*r  tern  in  Eq.  (3»l8).  For  Rx  , the 

-r 

close  cancellation  between  "curvature"  contributions  shows  that  the  7-r 

term  is  quantitatively  Important  for  Rj^  . On  the  other  hand,  our  analysis 

shows  that  it  is  negligible  for  R(j  . 

We  have  explained  the  temperature  dependence  of  the  Hall 

effect  within  our  model  in  terms  of  changing  Fermi  surface-Brillouin 

87 

zone  wall  interactions,  as  suggested  by  Stringer,  et.  al.  for  Cd. 

Hie  changing  c/a  ratio  (Table  3»l)  shifts  the  zone  planes  with  respect 
to  the  Fermi  sphere.  For  example,  7qqo  decreases  with  temperature. 
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increasing  the  relative  size  of  the  free  electron-like  "lens".  The 
doninant  temperature  dependence  in  our  model  derives  from  the  reduction 
in  band  gaps  . Quite  simply,  increasing  the  temperature  reduces  Fermi 
surface  distortion  through  the  band  gap  parameter  g,  driving  the  'rT»n 
constants  in  the  negative  direction.  These  effects  appear  adequate  to 
explain  the  observed  temperature  dependence  (Fig.  3.1),  although  we 
cannot  rule  out  completely  the  effect  of  temperature  dependent 
anisotropic  relaxation  times,  particularly  for  Rj^« 

C.  Thermoelectric  Power 

Our  calculated  values  of  vs.  u provide  a good  check  on 
the  model  and  method  of  computation  used  in  this  section.  Fig.  4.3  is 
in  qualitative  agreement  with  Fig.  C.3  from  Ref.  22;  the  lower  average 
values  obtained  here  being  due  to  the  simplified  model  ve  used.  However, 
since  we  are  interested  primarily  in  the  energy  derivatives  of  the 
relaxation  times,  and  not  on  their  actual  values,  this  disagreement 
does  not  appear  serious. 

*<rj> 

We  find , and  consequently  s?, , to  be  negative,  contributing 

de  4:1 

toward  a positive  thermopower.  As  in  the  noble  metals,  this  is  due 
primarily  to  the  fact  that  the  form  factor  has  its  first  node  at  a 
wavenumber  less  than  2ky,,  so  that  the  scattering  cross  section  (depending 
on  F ) becomes  a minimum  at  this  point  and  then  increases  again  beyond 
the  node.  This  results  in  an  increase  in  back  scattering  (large  values 
of  | k*-k  | ) probability  as  energy  increases;  such  processes  dominate 
the  resistivity  in  the  high  temperature  (umklapp)  regime.  As  might 
be  expected,  this  tem  proved  to  be  very  sensitive  to  the  exact  shape 
of  the  form  factor  used. 
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Our  model  (Fig.  4.1)  obviously  provides  only  a very  rough 
treatment  of  the  effects  of  Fermi  surface  distortion,  since  it  includes 
only  the  missing  area  correction.  As  the  Fermi  surface  expands  outward 
(i.e. , € increases)  against  the  slanting  Brillouin  zone  planes,  the 
amount  of  missing  area  is  increased,  and  the  integrals  X.^  correspondingly 
reduced.  The  missing  area  is  concentrated  nearer  the  basal  plane,  as 
reflected  in  the  anisotropy  in  the  parameter  s^  (Table  4.1).  An  independent 
calculation  of  only  the  s^  terms,  but  taking  full  Fermi  surface  distortion 
into  account  using  the  2-OIV  model  surface  of  APPENDIX  B,  yielded 
* -5.4  and  slc  = +.37.  These  values  are  in  rough  agreement  with 
those  in  Table  4.1,  the  differences  pointing  out  the  approximate  nature 
of  the  models  and  the  sensitivity  of  the  calculation.  Nevertheless, 
the  s1  factor  does  qualitatively  account  for  the  experimentally  observed 
anisotropy  in  Also,  the  suras  of  the  s^  and  s2  factors  result  in 

values  of  S and  S which  appear  consistent  with  experimental  values 

& c 

at  lower  temperatures  (Fig.  4.4). 
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APPENDIX  A 


DISCUSSION  OF  FORM  FACTOR 


The  pseudopotential  method  in  the  theory  of  simple  metals 

2k 

la  discussed  in  detail  by  Harrison.  His  formulation  is  based  on  the 

Orthogonallzed  Plane  Wave  (OFW)  method  in  which  the  basis  functions 

for  expansion  are  plane  waves  whose  projections  onto  the  core  states 

have  been  eliminated  to  insure  orthogonality  between  core  and  conduction 

states.  Lumping  the  core  projection  operators  into  the  potential 

provides  a repulsive  term,  cancelling  part  of  the  coulomb  attraction 

term.  The  resulting  pseudopotential  is  weak  (allowing  the  use  of  first 

order  perturbation  theory)  though  non-local;  i.e.,  its  matrix  elements 

between  states  k and  k in  general  depend  on  k and  k,  and  not  Just  on 

q ■ I?  - ?"|.  In  the  calculation  of  transport  properties,  however,  we 

are  interested  only  in  scattering  on  a constant  energy  shell,  and  the 

error  introduced  by  treating  the  effective  potential  as  local  is  small. 

The  solid  line  in  Fig.  A.l  shows  the  form  factor  for  Zu^ 

calculated _by  a semi-empirical  method.  The  scheme  involves  interpolation 

between  points  determined  by  Stark  and  Folicov^ 1 from  their  analysis 

of  low  temperature  Fermi  surface  measurements,  i.e.,  de  Haas -van  Alphen 

effect  data.  Calculation  of  the  core  repulsion  contribution  to  the 

pseudopotential  is  also  described  in  Ref.  57«  A recent  form  factor 

84 

calculation  by  Greenfield  and  Wiser  agrees  quite  well  with  the  values 

which  we  used.  In  numerical  computation,  we  used  the  analytic  expression 

sin(Bq/kJ  2 

F(q)  - A — + C sin  (Dq/kJ  , (A.1) 

(‘lAy) 
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vith  A * -2.^9  rydbergs,  B » 1.85,  C ■ .046  rydbergs,  and.  D * 2.03. 

The  analytic  fora  fit  the  semi-empirical  values  to  within  about  .01 

rydbergs.  Also  shown  in  Fig.  A.1  is  the  model  fora  factor  used  by 

22 

Chan  and  Huntington  ~ in  the  calculations  described  in  APPENDIX  C. 

This  is  given  by  the  functional  fora 

F(q)  - F0(l  - q2/2.9  k£)  , (A. 2) 

with  Fq  = -.34  rydbergs. 
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APPENDIX  B 
FERMI  SURFACE  MODEL 

The  Fermi  surface  of  Zn  is  distorted  by  three  sets  of 
Brlliouia  zone  planes  whose  RLV's  are  the  vectors  Q^,  Q2  and  of 
Eqs.  (4.17).  (We  switch  to  the  notation  K^,  reserving  ^ for  electron 
phonon  umklapp. ) The  amount  of  distortion  caused  by  each  plane  is 
determined  by  the  coupling  strength  parameter,  g = k^lF^VCpK  (see 
Eq. (B.l)  ).  The  relevant  g values  for  Zn  at  650  K sore  given  below. 

1*1  fry)  g 

\ 1.596  .023 

itg  1.740  . 001 

1.914  .047 

Table  B.l  : Coupling  Strength  Parameters 

Or 

The  tabulated  values  include  a correction  for  the  Debye-Waller  factor, 
which  effectively  reduces  the  bund  gaps  with  temperature.  The  table 
reveals  two  important  points.  First,  the  wavenumbers  of  the  Brillouin 
zone  planes  that  intersect  the  Fermi  sphere  fall  in  the  region  where 
the  form  factor  (Fig.  A.1)  is  relatively  small;  i.e. , the  region  where 
cancellation  between  Coulomb  attraction  and  core  repulsion  terms  in 
the  pseudopotential  is  close.  This  results  in  small  band  gaps  and  a 
NFE  band  structure.  Second,  the  Kg  RLV  falls  very  near  the  zero  of 
F(q).  We  will  therefore  ignore  the  distortion  caused  by  this  set  of 
planes,  which  results  in  considerable  simplification. 

Fig.  B.l(a)  shows  the  Jones  (or  second  Brillouin)  zone  for 
Zn,  and  Fig.  B.l(b)  depicts  the  Intersections  of  the  Fermi  sphere  with 


(a)  Jones  Zone 


(b)  Fermi  Sphere  Intersections 


Fig.  B. I!  Intersection  of  Fermi  Sphere  With  Jones  Zone 
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the  planes  determined  by  the  sets  of  vectors  1C^  and  K^.  We  treat  the 
effects  of  distortion  of  one  pair  of  planes  at  a time,  as  illustrated 
in  Fig.  B.2.  In  this  2-OIV  model,  the  equation 


€_ 

*F 


+ 2 ^ 1 


£(7-M)  ±[(7-M)2  + 


(B.l) 


gives  the  electron  energy  near  the  Kth  zone.  The  plus  sign  in  Eq.  B.l 
goes  with  u >7*  Differentiation  gives  the  electron  velocity  near  the 
Kth  zone: 


H? 


(B.2) 


Jb-V)2  + 62 


m 


(B.3) 


(B.4) 


Eq. (B.l)  also  reveals  that  for  the  slanting,  or  RLV's,  the  alues  of 
7 and  g are  such  that  the  region  corresponding  to  u>7  is  not  occupied. 
This  results  in  a considerable  amount  of  missing  free  area  from  the 
Fermi  surface,  as  illustrated  by  the  shaded  areas  in  Fig.  B.3.  Fig. 

B.3  also  illustrates  the  division  of  the  surface  into  regions  in  which 
the  effect  of  a single  pair  of  planes  is  considered.  In  the  computations, 
each  of  the  regions  bounded  by  dotted  lines  was  approximated  by  a region 
of  equal  area,  but  cylindrical!^  symmetric  about  the  appropriate  RLV. 
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SOLUTION  OF  THE  BOLTZMANN  EQUATION 

22 

We  present  here  a summary  of  the  Chan  and  Huntington 
solution  of  the  Boltzmann  equation  (Eq.  (4. 3)  ) for  the  vector  mean 
free  path  X.  Briefly,  they  write  each  component  of  T in  the  form 
X^  = , and  then  develop  an  iterative  scheme  for  the  relaxation 

time  T^.  The  transition  probability  W(k,k*)  is  dominated  by  the 
electron-phonon  interaction  in  the  high  temperature  regime.  Experimental 
phonon  dispersion  curves  are  shown  in  Fig.  C.l.  These  are  used  to 
derive  an  approximate  set  of  dispersion  curves,  each  applicable  to 
umklapp  processes  involving  a particular  RLV.  (Normal  processes  are 
handled  similarly. ) These  fitted  dispersion  curves  are  shown  in  Fig. 

C.2. 

Results  for  <n>  and  <T*>  (corresponding  to  the  "5  field 

C Ql 

being  applied  parallel  and  perpendicular  to  the  c-axis,  respectively) 
are  shown  in  Fig.  C.3.  The  phonon  dispersion  curves  show  that  acoustic 
velocities  are  lower  for  modes  polarized  perpendicular  to  the  basal 
plane.  As  a result,  the  electron  scattering  is  greater  for  the  umklapp 
processes  parallel  to  the  c-axis,  and  the  relaxation  time  is  thus 
shorter  for  conduction  in  the  c direction.  However,  distortion  of  the 
Fermi  sphere  and  the  resulting  missing  area  discussed  in  APPENDIX  B 
has  a greater  effect  on  conduction  in  the  a direction.  The  Fermi 
surface  distortion  compensates  quite  closely  the  elastic  anisotropy, 
giving  a nearly  isotropic  conductivity  in  good  agreement  with  experiment. 
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